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ABSTRACT 

Context. Orion-KL is a remarkable, nearby star-forming region where a recent explosive event has generated shocks that could have 
released complex molecules from the grain mantles. 

Aims. A comparison of the distribution of the different complex molecules will help in understanding their formation and constraining 
the chemical models. 

Methods. We used several data sets from the Plateau de Bure Interferometer to map the dimethyl ether emission with different arcsec 
spatial resolutions and different energy levels (from Eup=18 to 330 K) to compare with our previous methyl formate maps. 
Results. Our data show remarkable similarity between the dimethyl ether (CH3OCH3) and the methyl formate (HCOOCH3) distri- 
butions even on a small scale (1.8''x0.8'' or ~ 500 AU). This long suspected similarity, seen from both observational and theoretical 
arguments, is demonstrated with unprecedented confidence, with a correlation coefficient of maps ~ 0.8. 

Conclusions. A common precursor is the simplest explanation of our correlation. Comparisons with previous laboratory work and 
chemical models suggest the major role of grain surface chemistry and a recent release, probably with little processing, of mantle 
molecules by shocks. In this case the CH3O radical produced from methanol ice would be the common precursor (whereas ethanol, 
C2H5OH, is produced from the radical CH2OH). The alternative gas phase scheme, where protonated methanol CH30H2^ is the 
common precursor to produce methyl formate and dimethyl ether through reactions with HCOOH and CH3OH, is also compatible 
with our data. Our observations cannot yet definitely allow a choice between the different chemical processes, but the tight correlation 
between the distributions of HCOOCH3 and CH3OCH3 strongly contrasts with the different behavior we observe for the distributions 
of ethanol and formic acid. This provides a very significant constraint on models. 
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1. Introduction 

The Orion protocluster region is the closest high mass star for- 
mation region to the Sun. We adopt a distance of 414+7 pc 
( |Menten et al.||2007|, consist ent within the error bars with the 
value s of 389 pc of^Sandstrom et al. (2007 ), 437 pc of Hiro^ 
etaL] ( [2007j ) and 419 pc of Kim et al. (2008). This region is re- 
markable because of the presence of high speed shocks which 
appear to be generated in the center of the Kleiman-Low in- 
frared nebula Orion-KL. These are reminiscent of a past explo- 
sive event and are most cl early seen in the H2 2.12 jim emis- 
sion( Allen & Burton|1993 ). Erom VLA proper motion measure- 
ments of two strong centimetric radiosources, it was proposed 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



that a very unique phenomenon had taken place some 500-1000 
years ago: the close encounter, or collision, of two or more rather 
massive stars. The objects involved in such a dynamical interac- 
tion could have included t he Becklin-N eugebauer obj ect (BN) 
and sources I and n (Gomez et al.|2QQ5j [Rodrigue z et al.|2QQ5[ 
|Goddi et"ar||2Qll[ |Nissen et al.||2012| rTraces of this explosive 
patt ern have recent ly been observed in the CO distribution by 
Zapata^ar](|2009|. 



This very recent and energetic event at the heart of the nebula 
provides unique conditions for the study of a rich interstellar 
chemistry. Many molecules could have been released from the 
grain mantles because of dust heating or sputtering by multiple 
shocks, especially the relatively large molecules for which no 
gas-phase-only formation route is satisfactory. 

Erom a set of twelve Plateau de Bure observations, we 
have recently analyzed the deuterated methanol isotopologues 
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Table 1. Main parameters of the IRAM Plateau de Bure interferometer data sets where transitions of dimethyl ether are 
detected. 



Bandwidth 


HPBW 


Spectral resolution 


Flux conversion 


Synthesized beam 


Pixel size 


(GHz) 


n 


(MHz) (kms-i) 


(1 Jy bearn"^) 


c X PA n 


C X -) 


80.502 - 8O.574R 
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0.625 2.33 


4.6 K 


7.63 x 5.35 15 


0.80 X 0.80 


101.178- 101.71T] 


50 


0.625 1.85 


15.8 K 


3.79 X 1.99 22 


0.50 X 0.50 


203.411 -203.485^ 


25 


0.625 0.92 


7.0 K 


2.94x 1.44 27 


0.32 X 0.32 


223.408 - 223.941^ 


23 


0.625 0.84 


17.3 K 


1.79 x0.79 14 


0.25 X 0.25 



^ Observations centered on (0^72000 = 05^^35'" 14H6, 6j2ooo = -05°22'3a'59) with a velocity of 6 km s"^ with respect to the LSR. 
^ Observations centered on (ofysooo = 05^35"^ 14!20, 6j2ooo = -05°22'36r00) with a velocity of 8 km s"^ with respect to the LSR. 



CH2DOH and C H3QD (|Peng et al.|[2QT2l ) and the methyl for- 
mate HCOOCH3 ( |Favre et al.|2 011a) emission in Orion KL. The 
latter work suggests a rather close morphological relation in sev- 
eral places of the Orion KL nebula between excited H2 emission 
at 2.12 Jim and methyl formate. 

Even though both methyl formate and dimethyl ether are 
oxygen-bearing molecules of medium complexity and both can 
be produced from gaseous methanol (e.g. Peelers et al.||2006| 



Table 2. Detected and blended transitions of dimethyl ether 
observed with the Plateau de Bure Interferometer toward 
Orion-KL. 



and/or from grain surface chemistry (e.g. Garrod et al. |20Q8 ) 



these two species are spectroscopically different, and their spec- 
tra behave differently with physical temperatur e. Furthermore, 
since the ear ly days of dimethyl ether studies ( Sn yder etal. 1974[ 
Clark et al.||1979| ) there was a hint of limited departures from 
LTE for dimethyl ether. Recently, we have carried out observa- 
tions at the EVLA of the KK-uK+i) = 6(1, 5) - 6(0,6) EE tran- 
sition of dimethyl ether at 43.4475415 GHz ( |Favre et aL|2011b| ) 
which shows that the distribution is very similar to that of methyl 
formate. 

In this paper we present the analysis of several dimethyl ether 
lines present in our Plateau de Bure data set. The maps of the 
CH3OCH3 emission at different frequencies and energy levels 
are shown in Sect. [3] In Sect. |4] we compare the spectra at the 
five main emission positions to synthetic spectra based on the 
temperatures derived from our previous methyl formate analysis. 
Finally, we discuss in Sect. |5] the similarity in spatial structure 
and velocity between our high resolution maps of dimethyl ether 
and methyl formate and examine its implication on the different 
chemical models of formation of these species, in particular in 
comparison also with our ethanol and formic acid maps. 
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D: detected. B: blended. PB: partial blend. 



2. Observations 

The data set consists in twelve data cubes obtained with the 
IRAM Plateau de Bure Interferometer, the parameters of which 
are presented in Favre et al. (2011a; see their Table 1). We used 
the GILDAS packag^for data reduction. The continuum emis- 
sion was subtracted in the data cubes by selecting line-free chan- 
nels as judged by careful visual inspection, discarding any con- 
taminated channels. Finally, we cleaned the data cubes, channel 
by channel, using the B. G. Clark method ( |Clark|1980) ). 

Table[T]presents the parameters of the four cubes where tran- 
sitions of dimethyl ether are detected. Spatial resolution ranges 
from 1.79'' x 0.79'' to 7.63" x 5.35" and spectral resolution 
from 0.84 to 2.33 km s'^ 



3. Dimethyl ether (CH3OCH3) frequencies and maps 

3.1. Dimethyl ether frequencies 

The data set includes several dimethyl ether (DME) lines. Table 
[2] lists the DME transitions, present in our data and detected or 



blended, taken from the CDMS database p] ( jMuUer et all2001 
"2005 ) up to Eupper ^ 650 K and based on the recent work of 
Endre s et al. |( |2009| ). The non detected transitions correspond to 
lines too faint to be detected (Syu^>6.5 for the detected transitions 
and SyL/^<4.1 for non detected). The detected transitions cover 
the energy range 19 to 330 K. 

Dimethyl ether is an asymmetric top molecule with two 
equivalent methyl groups undergoing large amplitude motions 
along the CO-bond. The two internal rotations cause a splitting 



1 http://www.iram.fr/IRAMFR/GILDAS 
^ http : //www. astro .uni-koeln. de/cdms 
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Fig. 1. Dimethyl ether integrated intensity maps obtained with 
the Plateau de Bure Interferometer: these panels show the sum 
of the emission of the transitions at 80.54 GHz, 101.56 GHz, 
203.42 GHz and 223.41 GHz. The first contour and the step 
are 0.5 Jy/beam, 0.2 Jy/beam, 2 Jy/beam and 0.3 Jy/beam at 
80.54 GHz, 101.56 GHz, 203.42 GHz and 223.41 GHz respec- 
tively. Note that the 203.42 GHz transitions are blended with 
the H2^^0 line at 203.407 GHz around (a/2000 = 05^35^^14^.5, 
^/20oo = -05°22'34''). The BN object position is (a/2000 = 
05^35^ 14U 094, ^^2000 = -05°22' 22^724), the radio source I po- 
sition is (Qf/2000 = 05^35^^14^.5141, ^^2000 = -05°22'3a.'575), and 
the IR source n posi tion is (0^/2000 = 05 ^35^14^.3571, ^/2ooo = 
-05°22'32r719) from |Goddi et"aLl ( [2QTT] ). 



of each rotational level into four substates AA, EE, EA, and AE 
(see Sect. |4] for line strengths and frequency separation of the 
multiplets). 

3.2. The dimethyl ether emission maps 

The maps of the dimethyl ether molecule, CH3OCH3, allow us to 
trace the spatial distribution of one major oxygenated molecule 
in Orion-KL. Figure [T] shows the CH3OCH3 emission measured 
at diff'erent wavelengths; note the diff'erent angular resolutions of 
the various data sets (listed in Table [T]). The distribution shows 
the same extended, V-shaped molecular emission linking the ra- 
dio sour ce I to the BN obje ct as previously observed in methyl 
formate ( Favre et al.|2011a| ). The highest spatial resolution map 
at 223.41 GHz (1.79''x0.79'0 is presented in Fig. [2] We have 
marked the Hot Co re and Compact Ridge positions (see e.g. 
Beuther et al.||2005"] ). The position of the main emission peaks 
identified in methyl formate, labeled MFl to MF28, are indi- 
cated. Most of these correspond to the dimethyl ether peaks as 
well within the beam size. 

The distribution of CH3OCH3 appears very similar to the 
distribution of HCOOCH3 in the velocity range where the main 
component lies between 6 and 9 km s"^ and a north-south linear 
structure shows up at higher velocities (9-12 km s"^). 

The relative intensities of the diff'erent spatial emission peaks 
depend on the upper state energy of the transitions. Figure |3] 
shows the CH3OCH3 emission for transitions with diff'erent en- 
ergy levels ranging from Eup=18 to 330 K. The 43.47 GHz map 
is that obtained with the EVLA (Fa vre et al.||2011b ). As previ- 
ously observed for the methyl formate in |Favre et aLp 201 la ), the 
emission of dimethyl ether at the MF2 position ("Hot Core SW" 
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Fig. 2. Dimethyl ether integrated intensity maps obtained with 
the Plateau de Bure Interferometer (sum of emission at 
223.409 GHz and 223.412 GHz between 5 and 12 km s'^). The 
bottom image is a blowup of the Hot Core/Compact Ridge map 
area. The beam is 1.79'' x 0.79''; the level step and first con- 
tour are 3.8 K km s~^ Th e position of the millimeter source 
MM23 ( Eisner et all|2Q08l ) is also indicated . The main diff^er- 
ent HCOOCH3 emission peaks identified in F avre et aL]p01 la| ) 
are marked by a cross and labeled MFnumber- Note that the 
dimethyl ether distribution is very similar to that of methyl for- 
mate (cf. Fig. 4 in [Favre et al.|201 la| and see Sect. |5.1.l] ). 



in 'F riedel & Snyder|[2008| ) becomes stronger at higher energy 
levels. 

The comparison of our interferometric spectrum at 
101.56 GHz to the single-dish spectrum taken with the IRAM 
30m (J. Cernicharo, private comm.) shows that little flux is miss- 
ing due to filtering on a 3" scale (see Fig.|4]). 

4. Temperatures and column densities 

The four torsional forms AA, EE, AE and EA have weights of 
6:16:2:4 for ee, 00 levels and 10:16:4:6 for eo and oe levels 
( Turner! 1991 Lovas et aL]|1979| ). Most of the time the spectral 
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Table 3. Positions of the main HCOOCH3 emission peaks ob- 
served with the Plateau de Bure Interferometer toward Orion-KL 
as identified in Favre et al. ( 2Qllaj . These peaks also coincide 
with the main CH3OCH3 peaks. 
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Fig. 3. CH3OCH3 intensity maps integrated in velocity between 
5 and 12 km s"^ The line frequency and the upper state energy 
are indicated on each plo t. The CH3OCH3 4 3.47 GHz map was 
obtained with the EVLA (Favre et al.|20 11 b|). The dimethyl ether 
emission is stronger towards the compact ridge than the hot core 
SW position for low upper energy transitions, while the opposite 
is the case for high upper state energies. 



resolution of the radiotelescopes does not allow separation 
of the AE and EA transitions, so that a symmetrical triplet is 
observed. In that case and if the lines are optically thin, the 
observed dimethyl ether triplet (AA, EE, AE-hEA) should then 
be either 6:16:6 or 10:16:10 in relative intensity. For part of the 
triplet observed at the edge of the bandwidth at 223.4 GHz we 
find a ratio of 1 :2.7: 1 within the errors. This is in agreement with 
a low optical depth. However for the transitions at 80.53 GHz, 
the EE line in the middle is much too weak with respect to 
the other lines. This EE line is stronger in the spectra taken 
with single-dish radiotelescopes ( [Johansson et a l. 1985; Turner 
|1989| ), so that this could be caused by an observational rather 
than an excitation eff'ect. 

For the transitions at 101.56, 203.42 and 223.41 GHz, 
we have produced synthetic spectra assuming local thermo- 
dynamic equilibrium (LTE), using the line parameters derived 




1.0155 10 1.01555 10 1.0156 10 1.01565 10 1.0157 10 

Frequency (MHz) 

Fig. 4. Dimethyl ether spectrum at 101.56 GHz observed with 
the IRAM 30m radiotelescope (in black, J. Cernicharo, private 
comm.) overlaid on the PdBI spectrum convolved to the same 
spatial resolution (in red). 



from our methyl formate analysis toward the five main emission 
peaks whose positions are given in Table [3] We first used the 
HCOOCH3 velocity v, line width Avi/2 and temperature Trot to 
derive the NCH3OCH3 column density; then we have slightly ad- 
justed the parameters for a better fit of the dimethyl ether spec- 
tra. This was done with our own routines and XCLASS Fig. 
[5] shows synthetic spectra overlaid on the observed spectra and 
Table |4] lists the parameters for the angular resolutions of the 
101 GHz (3.63'' X 2.26'0 and 223 GHz (1.79'' x 0.79") ob- 
servations. The parameters used for the 203 GHz observations 
(angular resolution of 2.94" x 1.44") are intermediate values. 
The calculated opacities (estimated from the ratio of the bright- 
ness temperature to the rotational temperature) are generally less 
than 0.2 and at most 0.4 at MFl. 

We derived abundances relative to methanol using the col- 
umn densities of Peng et al. ( 2012 ). We find values between 0.05 
and 0.2. This result assumes that both gas components are co- 
extensive and in LTE at the same temperature, which might be 
questionable if the diff'erence in the derived temperature between 
dimethyl ether and methanol is significant. This question is re- 
lated to the dimethyl ether formation route and is discussed later 
in SectB 



5. Discussion 

5.1. Similarity of dimethyl ether and methyl formate 
HCOOCH3 maps 

Most maps of molecular emission in Orion KL obtained in pre- 
vious studies show a huge diversity in morphology. On the con- 
trary we show here the very striking similarity which appears 
between our high resolution maps of dimethyl ether and methyl 
formate ( Favre et al.|2011a| . The high degree of correlation we 
find reinforces the need to explain the similarity already noted 
between these species in previous studies, and probably points 
towards a similarity of their formation paths, with a common 
precursor for both species. 

First we present the correlation of dimethyl ether with 
methyl formate observed in this work and previous results for 
these molecules. We then briefly review the models of gas phase 



^ This research made use of the myXCLASS program 
(http s ://w w w. astro . uni-koeln . de/proj ects/schilke/XCL AS S ) , 
which accesses the CDMS (http://www.cdms.de) and JPL 
(http://spec.jpl.nasa.gov) molecular data bases. 
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Table 4. CH3OCH3 synthetic spectra parameters (velocity, line width, beam- averaged temperature and column 
density) which reproduce best the PdBI spectra at the emission peaks MFl to MF5 for the angular resolutions of 
the 101 GHz (3.63'' x 2.26'') and 223 GHz (1.79" x 0.79") observations. Two velocities components are clearly 
distinguished at MFl, MF4 and MF5. 
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101 GHz 



203 GHz 



223 GHz 





MF3- 




MF4- 



40 20 -20 
Velocity (km/s) 



10 -10 -20 
Velocity (km/s) 



MF5- 



10 
Velocity (km/s) 



Fig. 5. CH3OCH3 synthetic spectra (in red) overlaid on the ob- 
served spectra toward the MFl to MF5 positions. The 101 GHz, 
203 GHz and 223 GHz observations are in the left, middle and 
right columns respectively. The observed spectrum at 223 GHz 
is at the edge of the bandwidth. 



and/or ice mantle chemistry proposed to form these spe cies, con- 
centrating on the most recent works (Laas et al. 201 1 ; NeiUetal. 
Bisschop et al.|2007[|Qberg et al.,2Q09,,2Q10j and exam 
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ine the compatibility of various hypotheses with the present data, 
as well as with our ethanol and formic acid maps made from the 
present data set. 

5.1 .1 . PdBI maps and correlation diagrams 

When we compare the dimethyl ether map with the methyl for- 
mate map obtained with the same data set and from transitions 



with close energies (Fig. [6]), these show a striking overall sim- 
ilarity: 1) both molecules are present in the 6-11 km s"^ veloc- 
ity range and show no detectable emission at 5 km s"\ 2) both 
molecules show a general V-shaped distribution, 3) they both 
have their strongest peak at the Compact Ridge for transitions 
of low Eu energy range (MFl position), 4) they both have a sec- 
on d peak (MF2 position) a t the position named Hot Core SW 
by |Friedel & Snyder] pOQ8 ), a few arcseconds below the contin- 
uum Hot Core peak, and which is stronger than MFl for high 
Eu energy range transitions (see Sect. |3.2| ), 5) they both have a 
northern extension toward BN around 10 km s"\ 6) many other 
secondary peaks are similar. 

The correlation between the two maps is however not 100%. 
Faint emission is present for both species but differ somewhat in 
shape at other places or is displaced by 1-2''. One should note 
however the influence of the cleaning procedure necessary to re- 
construct the distribution. 

The similarity between dimethyl ether and methyl formate 
maps can be displayed in a more quantitative way. In Fig. |7] 
the intensity of each pixel of the dimethyl ether emission line 
(Ej^=330 K) at 223.41 GHz and methyl formate emission line 
(E^=305 K) at 223.534 GHz are compared (note that there is 
some correlation between pixels since the pixel size is about 3 
times smaller than the synthesized beam). The correlation coef- 
ficient of Bravais-Pearson is 0.8. By contrast we show the same 
type of diagram with formic acid HCOOH (10(2,9)-9(2,8) tran- 
sition at 223.915 GHz with Eu=12 K) compared to the methyl 
formate transition at 223.500 GHz (11(4,8)-10(3,7), E^=50 K) 
and ethanol C2H5OH (23(7,17)-23(6,17), v^=0-l transition at 
223.629 GHz with E^=346 K) compared to the methyl formate 
transition at 223.534 GHz (18(5,14)-17(5,13), E^-305 K), for 
which the correlation coefficients are only 0.32 and 0.19 respec- 
tively. For reference, we also show the correlation of 2 transi- 
tions of methyl formate with same energy (223.465 GHz and 
223.500 GHz with E^=50 K); the correlation coefficient is 0.82 
- the methyl formate/dimethyl ether result is almost as high. The 
correlation of 2 transitions of methyl formate with diff'erent ener- 
gies (223.465 GHz and 223.534 GHz with E„=50 K and 305 K, 
respectively) is not as good (correlation coefficient of 0.76) due 
to differenc es in excitation which are apparent in Fig. [3] and in 
Favre et al. (2011a, see their Fig. 15). The estimated opacity for 
all the above transitions is less than 0.15, so that there is no no- 
ticeable optical depth effect. 
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This similarity of the dimethyl ether map with the methyl 
formate map is not observed with other species (see e.g. Guelin 
let al.,2QQ8, and our Fig.[8]for formic acid HCOOH). 



-5°22'20" - 



-5°22'40" - 
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Fig. 6. Comparison of the methyl formate map (left) at 
223.534 GHz (E^=305 K) and the dimethyl ether map (right) 
at 223.41 GHz (Ef^=330 K). The transitions are from the same 
data cube and the 1.79'' x 0.79'' beam is plotted in the bottom 
left corner. 




Fig. 7. A plot of the intensity of pixels in the methyl formate 
maps versus the intensity of the same pixels in the dimethyl 
ether (top left), ethanol (top right), formic acid (bottom left) 
and methyl formate (bottom right) maps. Care has been taken 
to use transitions with similar Eu energies. The methyl formate 
(Ef^=305 K) at 223.534 GHz is compared to the dimethyl ether 
transition (Ej^=330 K) at 223.41 GHz and to the ethanol transi- 
tion (E^=346 K) at at 223.629 GHz. The methyl formate transi- 
tion (Eu=50 K) at 223.500 GHz is compared to the formic acid 
transition (E^=72 K) at 223.915 GHz and to the methyl formate 
transition (E^=50 K) at 223.465 GHz. Dimethyl ether shows 
clearly the highest spatial correlation with methyl formate. 



5.1 .2. Dimethyl ether and methyl formate behavior in 
previous publications 

The similarity in the distribution between the two species has 
also been seen in single dish studies of Orion and of other re- 
gions. In Table [5] we selected only some of the numerous Orion 



KL line surveys based on two factors: 1) early papers from 
which the parallel between some species so called "compact 
ridge species" , and in some cases already the parallel dimethyl 
ether/methyl formate distribution, was progressively recognized 
with increasing evidence, and 2) papers providing the best line 
profile parameters, temperatures and column densities for both 
dimethyl ether and methyl formate. 

When we convolve our PdBI dimethyl ether data to a typical 
single-dish field of view of 30", the spectra can be relatively 
well fitted with one component of temperature T^t ~ 90 K and 
column density NCH3OCH3 



0.6x lO^^cm-^. 



Using Mx = 4.64 x 10"^^ x yu^ x 



to convert the column density of molecule X (averaged over a 
beam of diameter 6) into a mass (in solar mass M©), we find 
for dimethyl ether (jix = 46): Mdme ~ 8.6 x 10"^ M©, or about 

2.8 MEarth- 



To better compare our results with the single dish studies, 
we have computed the dimethyl ether mass in each case, and 
we display the value, R, relative to our result, in column 9 of 
Table [5] All values are withi n a factor of two w ith respect to 
ours except in the case of the |Lee & Cho| ( [2002] ) observations. 
No trend is seen with respect to the beam size, despite the large 
range of sizes; this implies that at least 50%, and perhaps almost 
all of the dimethyl ether seen in the ODIN wide beam (2.1') is in 
fact contained in the inner 30" which we have observed with the 
PdBI. There is no noticeable eff'ect either related to the range of 
upper level E^^ used in the various studies. 



The parallel between both molecules is especially striking 
on small scales in the present PdBI data. Previous interferomet- 
ric papers showed some hints of this tendency. jMinh et aL ( 1993 1 
present the first HCOOCH3 and CH3OCH3 maps, made with the 
NMA (5'.' 2 X 4'.' 2 resolution); they noted that the diff'erence be- 
tween the emission distributions of the two molecules may be 
due to the 10 K differenc e in upper energy of the observed tran- 
sition levels. Blake et al.| ( p^996) with OVRO at 1.3 mm (1'.'5 x 
I'.'O synthesized beam, and 1 MHz =1.3 km/s spectral resolu- 
tion) observed 18 species including methanol, methyl formate 
and dimethyl ether. No maps are shown for dimethyl ether, but 
there is a general trend to behave somewhat like HCO OCH3, 
CH3OH, H2CS and H2CO. IBeuther et all ( [20Q5| [2Q06| show 
many maps obtained with the SMA at 0.86 mm and 0.44 mm 
of 0-bearing species including dimethyl ether, methanol and 
methyl formate, but at diff'erent excitation levels and there is 
only a partial similarity of the maps. In their CARMA ob ser- 
vations at 1.3 mm (2'.'5 x a.'85), [Friedel & Snyder] ( [2008] ) ob- 
served both dimethyl ether (transition 3XEu=^l K) and methyl 
formate (transition at E^=120 K). The maps of the two species 
did not match perfectly, one peaking at IRc5 and the other at 
IRc6, but the CH3OCH3 line is blended and a blend with an N- 
bearing molecule would explain the observed strong emission 
at IRc6 and the Hot Core SW. Note also that no methyl for- 
mate is attributed to the Compact Ridge (while it is the strongest 
methyl formate peak in Favre et al.,2011aj ); but these authors de- 
fine the Compact Ridg e position d i ff'erent ly. The correlation is 
very clear in Fig. 3 of |Neill et al.| (|2011) from their CARMA 
observations; the 6'.'1 x 5'.'0 resolution is however insufficient to 
clearly separate the emission from the many individual methyl 
formate peaks. 
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Table 5. Dimethyl ether and methyl formate parameters from previous single-dish studies. 
Authors 



( 2009|) 



Telescope 


Beam 


Frequency 


Dimethyl ether 


Methyl formate 


R^ 








Trot 


N 


Trot 


N 


1 1 




n 


(GHz) 


(K) 


(lO^^cm-^) 


(K) 


(lO^^cm-^) 




GBT 


16 


43 


88 


5.8 






1.07 


Onsala 20m 


47 


72-90 


75 


0.25 


25-90 


0.05-0.25 


1.02 


NRAO 11m 


83 


75-115 


91 


0.13 


62 


0.66 


1.66 


TRAO 14m 


46 


138-164 


130-360 


0.5-3.6 


55-120 


0.6-4.4 


2-14 


NRAO 12m 


43 


130-170 


75 


0.6 






2.05 


FCRAO 14m 


40 


150-160 


109 


0.3 


52 


0.077 


0.89 


OVRO 10.4m 


30 


215-263 


63+5 


0.3 


90+10 


0.26 


0.50 


JCMT 


14 


334-343 


75 


1.49 


105 


2.45 


0.54 


CSO 


20 


325-360 


89+5 


1.8 


98+3 


1.5 


1.33 


CSO 


11 


607-725 


360 


3 


316 


13 


0.67 


CSO 


11 


795-903 


160 


2 






0.45 


Odin satellite 


2.r 


486-492/541-577 


112 


13 






0.73 


PdBI 


30 


101/203/223 


90 


0.6 






1 



Comment 







Goddi et al. 



Johanssonetl 
Turner (199 f{ 
Eee & Cho ( 2002]) 

Remijan et al. (2 0Q8| 

Ziurys & McGon agTe| ( |1993| ) 
makeet al. (1987j~ 
Sutton et al. (199^ 
Schilkeet al. (19^1 
Schilkeet al. ( ZOOTi 
Comito et al.|( [20n5 
Fersson et aq pOOV I 
TEis paper ^ 



5.1.2 



^ Normalized number of dimethyl ether molecules, see Sect. 
^ Assuming a 10'' source size. 
^ a second cooler (Trot= 22 K) component is seen in MF only with a lower column density (2.7 10^"^). 
^ Compact Ridge position. 
^ Adopted source size 5-6''. 

^ Our PdBI data convolved to a field of view of 30". 



5.1 .3. The abundance correlation of dimethyl ether with 
methyl formate 

Another important input is the value of the relative dimethyl 
ether/methyl formate abundances. Table [6] lists the abundances 
derived at the main positions in our maps for dimethyl ether 
with respect to methyl formate. Dimethyl ether is about 3-4 
times more abundant than methyl formate towards MF3, MF4 
and MF5. 

Several hypotheses can be considered to explain this corre- 
lation: 



- Dimethyl ether and methyl formate would emit at the same 
location because of a similar excitation of the transitions, 
i.e. a similar behavior with respect to temperature, density 
or IR radiation field. But there is no special similarity in the 
level structure of dimethyl ether and methyl formate, espe- 
cially as compared to other complex organics like ethanol. 
Furthermore, the physical parameters are not uniform over 
Orion KL: the temperature varies from 80 to 160 K for 
methyl formate, which aflTects the dimethyl ether/methyl for- 



mate line ratio. The critical densities are 

l4-7 ™-3 



10^ 



cm 



^and 

,4 



~ 10^"^ cm"-^ respectively, while nH2 varies from 10^ to 
10^ cm"^. And the IR absorption bands are dififerent, only 
the CH3-O bond is common. 

- These molecules are just undergoing sublimation at the same 
temperature. This is likely, but similar abundances in the ices 
are also required in that case. 

- The similar column densities reflect only the general struc- 
ture of Orion KL. However dust and other molecules have 
diflTerent patterns. 



This leaves a correlation of the molecular abundances in the 
gas phase as the most simple and likely explanation (see further 
discussion in Sect. 15. 2b. 



Table 6. CH3OCH3 abundance relative to HCOOCH3 at the 
emission peaks MFl to MF5 for the angular resolutions of the 
101 GHz (3.63'' X 2.26'0 and 223 GHz (1.79'' x 0.79") ob- 



servations. We used the Nhcooch3 values listed in Favre et al. 
( |2011a| ). 



Emission peaks 


NCH3OCH3/ 
(1.79" X 0.79") 


NHC00CH3 
(3.63" X 2.26") 


MFl ^ 
MFl ^ 


2.15 
0.22 


9 
4 


MF2"^ 


1.6 


3.5 


MF3 


3.3 


3.6 


MF4 




4 


MF5 




3.6 



Component at 7.5 km s ^ . 
Component at 9.2 km s~^ . 



5.1 .4. Dimethyl ether and methyl formate observations in 
other regions 

Single dish observations of several other sources have also noted 
the parallel between these two species. 

In their study of 7 mass i ve hot core regions spread across 
the Galaxy, I Bisschop et al.| p007| ) study the statistical corre- 
lation between the abundances of 13 organic molecules look- 
ing for evidence of grain surface chemistry. Methyl formate and 
dimethyl ether are mostly seen where T^ust is > 100 K. A strong 
correlation is found between the abundances of H2CO, CH3OH, 
C2H5OH, HCOOCH3 and CH3OCH3 and the correlation is even 
> 0.9 among the last three. 

Comparing the results of [Bisschop et al. ( |2QQ7 ) with our 
findings, we note that 1) if dimethyl ether appears to be strongly 
correlated with methyl formate, the correlation they observed is 
similarly high with ethanol, 2) Trot has a large scatter among 
sources of dimethyl ether, and in some sources may be lower or 
higher than for methyl formate. Beyond possible observational 
uncertainties (noise and confusion), one should note the diflTer- 
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ent scales of both studies: while we sample down to {1.19" x 
0.79'' or -500 AU) in the nearby Orion region with the PdB in- 
terferometer, the single dish results of Bisschop et al.|p0d7 ) are 
average values for the whole region at Tdust > 100 K which is 
1500-5300 AU in radius for their sources - the latter would corre- 
spond to a study of Orion over ~ 20'' in diameter, i.e most of the 
KL region. This would explain why ethanol is not distinguished 
from dimethyl ether and methyl formate, and the less tight cor- 
relation between dimethyl ether and methyl formate could be 
due to additional formation pathways/excitation conditions for 
a fraction of these species on a larger scale. In other words, the 
process of formation of methyl formate and dimethyl ether we 
observe in our case, probably linked to shocks, might not appear 
as "pure" and unique on a larger scale. 

It is interesting to note that in a recent multi- species in- 
terferometric study o f the young binary protostar IRAS 16293, 
J0rgensen et aLlpUTl ) have mapped CH3OCH3 and HCOOCH3. 



The two molecules appear to be similar, showing two small con- 
centrations of similar intensity around the A and B source com- 
ponents. 

An extensive study of -bearing species was performed in 
the Galactic Center clouds by |Requena-Torres et al.| ( |2006||2008| ) 
which led them to conclude to a "universal" mantle compo- 
sition in the v arious GC clouds; this com position is different 
in hot corinos ([Reiquena-Torres et al.| [2Q06) and in Hot Cores 
( Requena-Torres et al. |2008| ). They relate the gas phase abun- 
dances of these complex organic molecules to sputtering or 
evaporation due to shocks, a situation similar to some extent 
to Orion KL. In the 2006 study, 7 lines of methyl formate and 
2 of dimethyl ether were observed, with E^/k in the same 20- 
40 K range; dimethyl ether was not observed in their 2008 study. 
Note that the rotational temperatures derived in the 2008 study 
are rather cold (<15 K). Due to the large distance of the GC 
clouds the sources they observe are likely more extended than 
what we observe in Orion. These authors do not put any empha- 
sis on the dimethyl ether/methyl formate comparison, although 
the two molecules are noted to have similar abundances with re- 
spect to methanol. The grain surface scheme proposed (Fig. 9 in 
[Requena-Torres et al. 2 008 ) does not suggest any close relation 
between methyl formate and dimethyl ether. 

5.2. Formation of dimethyl ether and methyl formate 
5.2.1. Chemical models 



Since detailed presentations can be found elsewhere (e.g. Peelers 
leFal. 2006; Neill et al.|20TT||Herbst & van Dishoeck|2009| ) we 
summarize here only briefly diff'erent models. 

A first model of pure gas-phase ion-molecule chemistry was 
proposed (e.g. |Blake|198l) ) where the protonated methanol ion, 
CH3OH2+, produced both CH3OCH3 and HCOOCH3, through 
reactions with CH3OH and H2CO respectively, which natu- 
rally explained the suspected intimate chemical link between 
these species. However, a first difficulty for an interstellar chem- 
istry model is to produce a sufficient amount of these complex 
species. The pure ion-molecule chemistry proved unable to pro- 
duce the observed abundances. This led to suggest a role of 
surface chemistry and of the release of mantle molecules into 
the gas phase by some process. The models considered first 
injection of water, but injection of methanol itself, processed 
once in the gas phase by ion-molecule chemistry appeared to 
be required. A key species in this scheme remained protonated 
methanol CH30H2^ , and a close correla tion of methyl formate 
and dimethyl ether was predicted (e.g. Charnley et al. 1995[ 



Caselli et al. |1993 ). Most o f the more complex sp ecies pre- 
dicted to be abundant by the 'Charnley et al.| ( [1995| ) model re- 
main undetected, although one of these (not expected to be the 
most abundant one), ethyl formate HCOOC2H5, has recently 
been detected in Sgr B2 ( Belloche et al.|(2009j . This ice man- 
tle methanol release/gas-phase post-processing model could not 



produce enough methyl formate when calculations by |Horn et al. 
( 2004) showed that the reaction CH30H2^ + H2CO was not pos- 
sible at interstellar temperatures. 

Presently there are two diff'erent proposals: 1) new gas-phase 
ion-molecule formation routes where formaldehyde is replaced 
by formic acid (e.g. [Neill et al.| 2011); 2) a direct formation of 
methyl formate and dimethyl ether on the grains (e.g. |Garrod| 
let al.|2008l[Laas et al.|201l] ). 

5.2.2. How do these models confront our Orion dimethyl 
ether and methyl formate data 

In our view, the simplest and most natural explanation of the 
dimethyl ether-methyl formate correlation is the production 
of both species by a common precursor reacting with an 
abundant simple species such as either CH3OH, H2CO or 
HCOOH. However, we suggest another - less likely in our view 
- possibility: the convergence of the dimethyl ether and methyl 
formate abundances resulting from the long time evolution of 
many chemical reactions. 



The common precursor hypothesis: CHO radical (solid phase) 
vs CH3OH and CH3OH2+ (gas phase) 

If solid phase chemistry dominates, the common precursor 
would be the CH3O radical, as detailed in Oberg et al (2010). In 
the scheme of UV irradiated ices, photodissociation of CH3OH 
leads to two radicals CH3O and CH20H|^ The former is the pre- 
cursor of methyl formate and dimethyl ether which are formed 
respectively by reaction with HCO and CH3 on the grain. The 
other radical produced from methanol, CH2OH, leads to ethanol 
C2H5OH, ethylene glycol CH2OHCH2OH and glycolaldehyde 
CH2OHCHO. A strong prediction of this model would thus be 
that the 3 latter molecules should show a somewhat similar dis- 
tribution, which is likely to be diff'erent from that of methyl 
formate and dimethyl ether. Diff'erent conditions, either in the 
pre-explosion phase (a preencounter role of BN, distance to the 
Trapezium, diff'erent heating, e.g. due to the proximity of source 



^ The mantle chemistry scheme needs UV radiation to produce rad- 
icals. Is it possible i n Orion KL? Acc ording to our previous study on 
deuterated methanol ( |Peng et al.|2012) ), the composition of ice mantles 
might have been largely determined in an earlier phase of the cloud, 
well before the explosive event. During this earlier phase, the matter 
was conceivably more uniform and less dense, hence less opa que to UV. 
The present very high opacity of the cloud (up to A^ > 1000, Favre etal.| 
2011a) suggests that processing by external UV sources is now negli- 
gible except on a very thin external layer. I nside a dense cloud UV can 
also be produced by dissociative shocks (e.g.|Neufeld & Dalgarno|1989]) 
or by secondary electrons (e.g. [Gredel et al.|1989[ [Prasad & Tarafdar| 
|1983 |. The presence of at least one B star inside Orion KL (the BN ob- 
ject) could also be an internal UV source but which is unlikely to reach 
the compact ridge due to absorption. However, as we have no knowl- 
edge of the precise 3D structure of the cloud, we cannot completely 
exclude that there are holes in the gas and we ignore the 3D dust dis- 
tribution within Orion-KL. But the simplest origin for the UV mantle 
irradiation remains that it occurred at an earlier epoch. Note that radicals 
can be for med also in mantles also by irradiation by charged particles 
(see e.g. |Bennett & Kaiser|2007}|Herbst & van Dishoeck|2009f . 
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I), or after the explosive event (direct effect of the shock, or sub- 
sequent exposure to the photons of BN, even higher proximity 
of IR luminous source I) might explain the different importance 
of the CH2OH and CH3O molecules in the Hot Core and in the 
Compact Ridge. Indeed, our interferometric observations show a 
similar distribution for C2H5OH and CH2Q HCH2OH but differ- 
ent from that of HCOOCH3 and CH3OCH3 dGuelin etal.|2008l 
Brouillet et al. in prep.); however, CH2OHCHO is not detected 
( [Favre et al.|2011a| ). Note that some conversio n mechanisms be- 
tween CH2OH and CH3O are suggested by Cernicharo et al. 
( |2Q12| ) in their CH3O detection paper. 

If the molecules would be mainly produced in the gas phase, 
methyl formate could be produced from formic acid as sug- 
gested by |Neill et al.| ( 2QTT] ): the anti-correlation observed be- 
tween HCOOH and HCOOCH3 across most of the Orion-KL 
region (see Fig. [8]) is consistent with their model of recent gas- 
phase conversion. Two ion-molecule reactions involving the re- 
action of methanol and formic acid, where one of the reactants 
is protonated, could be viable interstellar reaction routes to form 
methyl formate. Either the methanol, CH3OH, or the protonated 
methanol, CH30H2^, would then be the common precursor to 
methyl formate and dimethyl ether. Indeed, the reaction between 
protonated methanol and neutral methanol produces protonated 
dimethyl ether and is considered as an important contributor to 
interstellar dimethyl ether formation. If methyl formate is effi- 
ciently produced in a region by the reaction between HCOOH 
and CH30H2^, dimethyl ether is produced as well in the same 
spatial region, since neutral methanol is also present. A pre- 
diction of these reactions is the presence of the less stable t- 
HCOOCH3, which is not yet detected in Orion, and the anti- 
correlation with HCOOH, consumed by the reaction. Another 
interesting test would be to detect directly CH30H2^, once its 
millimetric spectrum is known. 



-5°22'30' 



-5°22'40" 




14^ 



Fig. 8. Methyl formate channel map (black contours) at 8 km s"^ 
(sum of the transitions at 223.465 and 223.500 GHz, E^=50 K) 
overlaid on the formic acid channel map (colors) at the same 
velocity for the 223.915 GHz transition (E^=72 K). All data are 
taken from the same data set (1.79'' x 0.79'' resolution). 



Methanol maps (see maps in 'Peng et aL]|2012[ from the 
same data set) are somewhat different. Despite the fact methyl 



formate and dimethyl ether are chemically related to CH3OH, 
their distributions are not expected to correlate as tightly with 
methanol as with each other since: i) in the gas phase production 
scheme, the methyl formate/dimethyl ether precursor CH30H2^ 
can be produced with an efficiency varying across Orion-KL. In 
addition CH3OH, as the parent species, is partly consumed to 
some extent (5-20%) to form methyl formate and dimethyl ether 
(and could have other reactions); ii) in the ice mantle production 
scheme, some layering of the mantle is expected (see e.g. |HerbsT 
& van Dishoeck 2009), and only a fraction of solid methanol 
closest to the surface can be processed by UV. 

We suggest, however, another possibility which could be 
investigated for the gas-phase chemistry. Because the main 
problem to form HCOOCH3 from CH2OH2+ and H2CO is an 
energy barrier, could this obstacle be removed if one considers 
that the ions are strongly accelerated with respect to neutrals at 
some places of an MHD shock ? A few km/s would represent 
enough kinetic energy to overcome the reaction barrier, which 
is 128 kJ/mol (15 000 K or 1.2 eV) according to [Horn et aL 
(!2004|). New quantum calculations and/or laboratory work are 
required to test this hypothesis. 



An alternative hypothesis could be: convergence of methyl 
formate and dimethyl ether abundances with time. 
We cannot exclude the convergence with time of the chemical 
network to a same methyl formate/dimethyl ether ratio value 
for a large range of initial conditions and parameters. (Note that 
such a convergence is not always present; for instance chemical 
systems may oscillate.) No common precursor is needed in this 
case. The convergence should be achieved in a time short enough 
(at least < 1-10 Myr , the probable maximum age of molecular 
clouds in OMC 1). This hypothesis is less likely in our view (cf 
the probable consequences of recent shocks). 

6. Conclusions 

We have studied the distribution of the complex 0-bearing 
molecule dimethyl ether CH3OCH3 at medium and high angular 
resolution (6.5'' - V) using various sets of interferometric data 
from the IRAM Plateau de Bure Interferometer. Our main results 
and conclusions are the following: 

1. This data set includes 4 well detected lines from £^^=19 K to 
Ej^=330 K and other partially blended or blended lines. 

2. The most intense emission arises from the Compact Ridge at 
the methyl formate peak position (MFl, Favre et al.'201 1a|) ; 
the second most intense a t the Hot Core SW p osition (Friedel| 
|& Snyder|2008| or MF2 ( [Favre et al.|20TTa| . 

3. We have used the temperature derived from our previ- 
ous methyl formate study for the five main peaks and we 
have deduced CH3OCH3 column densities assuming LTE. 
Temperatures cover the range 80 to 170 K, and column den- 
sities 1.5x10^6 to 3.5x10^'^ cm-^ 

4. The abundance of CH3OCH3 relative to methanol, CH3OH, 
is in the range 0.05 to 0.2 depending on the positions and 
assuming identical spatial distribution and temperature for 
both molecules. When observed on the same spatial scale of 
3.6'' X 2.3'', the abundance of CH3OCH3 relative to methyl 
formate, HCOOCH3, is in the range 3.5 to 9. 

5. We observe a very good correlation of the spatial distribution 
of the methyl formate and dimethyl ether emission. We show 
that it is most likely due to a correlation of their abundances. 
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6. The dimethyl ether emission follows the 2. 12 jim H2 distribu- 
tion, as does the methyl formate distribution. Shocks seem to 
be related to the presence of these species, possibly because 
of the release of molecules from grain mantles. 

7. A common precursor to dimethyl ether and methyl formate 
appears the simplest explanation to the observed correlation. 
In the gas phase, the precursor CH30H2^ would re act re- 
spectively with CH3OH and, as suggested by Neill et al.| 
( 2Q11| ), with HCOOH. The reaction of protonated methanol 
CH30H2^ with H2CO and CH3OH is excluded by quantum 
calculations in normal interstellar conditions. We speculate 
whether the reaction of high speed CH30H2^ ions (a few 
km/s to > 20 km/s) might change this conclusion. 

On the other hand, if methyl formate and dimethyl ether are 
already produced in the ice mantle, CH3OH and/or methoxy 
radical CH3O seem to be the common precursor. 

8. We observe a different distribution of ethanol C2H5OH (and 
ethylene glycol CH2OHCH2OH, in prep.). In the grain man- 
tle chemistry scheme this could result from different produc- 
tion rates or schemes of CH2OH vs CH3O radicals across 
Orion-KL. 

Furthermore we confirm the anti-correlation of methyl for- 
mate with formic acid, HCOOH, found by |Neill et al.] ( [20TT] ). 
This in favor of their model of production of methyl formate 
from gas phase HCOOH and CH3OH2+ where CH3OH2+ is 
formed from methanol released from grains. 

High resolution mapping brings new insight to the compari- 
son of complex organic molecules and the understanding of their 
formation. The HCOOCH3 and CH3OCH3 tight correlation, and 
the different behavior of ethanol we observe are a strong con- 
straint for future chemical models. 

We identify the need for the following key studies to con- 
strain further the chemistry of dimethyl ether and 0-bearing 
complex species: search for mantle species produced by the 
CH2OH radical; search for CH3OH2+; maps of t-HC00CH3 and 
HCOOH; theoretical/laboratory studies of protonated methanol 
reaction at high speed (shocked) on H2CO and CH3OH. 
The conditions leading to the convergence of the dimethyl 
ether/methyl formate ratio in a complex chemical network could 
also be investigated. 

Shock related molecule release from grains and shock in- 
duced chemistry of complex molecules need to be modeled. The 
issue of the efficiency or lack of such chemical conditions to 
transform molecules in a short time lapse (<10^ yr) is criti- 
cal in the discussion of the merits of the two methyl formate 
and dimethyl ether production schemes (gas phase production 
vs simple release); this will be addressed in a future work. 
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